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ABSTRACT 



We report on the detection of a very rapid TeV gamma-ray flare from 
BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging 
Telescope Array System (VERITAS). The flaring activity was observed dur- 
ing a 34.6-minute exposure, when the integral flux above 200 GeV reached 
(3.4 ± 0.6) x 10~ 6 photons m~ 2 s _1 , roughly 125% of the Crab Nebula flux mea- 
sured by VERITAS. The light curve indicates that the observations missed the 
rising phase of the flare but covered a significant portion of the decaying phase. 
The exponential decay time was determined to be 13±4 minutes, making it one of 
the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spec- 
trum of BL Lacertae during the flare was soft, with a photon index of 3.6 ± 0.4, 
which is in agreement with the measurement made previously by MAGIC in a 
lower flaring state. Contemporaneous radio observations of the source with the 
Very Long Baseline Array (VLBA) revealed the emergence of a new, superluminal 
component from the core around the time of the TeV gamma-ray flare, accom- 
panied by changes in the optical polarization angle. Changes in flux also appear 
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to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of 
the flare, although they are difficult to quantify precisely due to sparse coverage. 
A strong flare was seen at radio wavelengths roughly four months later, which 
might be related to the gamma-ray flaring activities. We discuss the implications 
of these multiwavelength results. 

Subject headings: active galactic nuclei: individual (BL Lacertae) - gamma-rays: 
individual (VER J2202+422) 



Introduction 



Blazars form a subclass of active galactic nuclei (AGN) that feature a relativistic jet 
pointing roughly at the observer. They are known for being highly variable at all wavelengths. 
In the most extreme cases, the timescales of gamma-ray variability can be as short as a few 
minutes at very high energies (> 100 GeV; VHE). Such va riability has been d etected in 
several BL Lacerta e objects (BL Lacs), including Mrk 421 (IGaidos et al.l 119961 1. Mrk 501 



([Albert et al.l 120071 ). and PKS 2155-304 (lAharonian et al.l 12007), and more recently in the 



flat-spectrum radio quasar (FSRQ) PKS 1222+21 (lAleksic et al.ll201ll ). The rapid variability 
poses serious challenges to the theoretical understanding of gamma-ray production in blazars. 
On the one hand, rapid gamma-ray variability implies very compact emitting regions that 
can be most naturally associated with the immediate vicinity of the central supermassive 
black hole. On the other hand, the regions must be sufficiently outside the broad-line regions 
(BLRs) that gamma rays can escape attenuation due to external radiation fields (which, for 



(Ghisellini & Tavecchio 


20081 


Giannios et al. 


2009 


Nalewajko et al. 


20121: 


Naravan & Piran 


|2012|). 



2009: Tavecchio et al. 2011; Barkov et al. 2012 



The spectral energy distributions (SED) of blazars show two characteristic peaks, with 
one in the infrared (IR) - X-ray frequency range and the other in the MeV - TeV gamma- 
ray range, respectively. The lower-energy peak is believed to be associated with synchrotron 
radiation from relativistic electrons in the jet, and the higher-energy peak with inverse- 
Compton radiation from the same electrons in leptonic models; the situation is more complex 
in hadronic models. Going from high-power quasars to low-power BL Lacs, the peaks shift 
systematically to higher frequencies. Most of the known TeV gamma-ray blazars are BL 
Lacs. They have been historically divid ed into high-frequency-peaked BL Lacs (HBLs ) and 
low- frequency-peaked BL Lacs (LBLs) (jPadovani fc Giommilll995l : iFossati et al.l 119981 1 . BL 
Lacertae, the archetypical source of the class, is an LBL in this classification scheme. 



-5 - 



B L Lacertae (also k nown as 1ES 2200+420) is an AGN located at a redshift of z = 
0.069 ( IMiller et al.lll978[ ). In 1998, the Crimean Obse rvatory reported a d etection of the 
source at >100% of the Crab Nebula flux above 1 TeV (INeshpor et al.ll200ll ). Subsequently, 
the MAGIC Collaboration reported another detection during an activ e state in 2005, bu t 



at a much lower flux level (only about 3% of the Crab Nebula fl ux) (lAlbert et al. 



2007). 



Triggered by activities seen with the Fermi LAT ( Cutini et aU2011) and AGILE (Piano et al. 



20111) at GeV gamma -ray energie s, as well as in the op tical (jLarionov et al.ll201ll ). near-IR 



(jCarrasco et al.ll201ll ). and radio ( lAngelakis et al.ll201ll ) in 2011 May, we began to monitor 
BL Lacertae more regularly at TeV gamma-ray energies with VERITAS. In this work, we 
report the detection of a rapid, intense VHE gamma-ray flare from the direction of the source 
on MJD 55740 (2011 June 28), as well as the results from the multiwavelength observations 
that were conducted around the time of this flare. 



2. Observations and Data Analysis 

2.1. Very High Energy Gamma Ray 

VERITAS is an array of four 12-meter imaging atmospheric Cherenkov telescopes lo- 
cated in southern Arizona. Each telescope is equipp ed with a focal-pla ne camera with 499 



photomultiplier tubes, covering a 3.5° field-of-view ( IHolder et al.l 120081 ) . VERITAS is sen- 
sitive to VHE radiation in the energy range from ~100 GeV to ~30 TeV, being capable 
of making a detection at a statistical significance of 5 standard deviations (5 a) of a point 
source of 1% of the Crab Nebula flux in ~25 hours. 

Prior to the intensified monitoring campaign with VERITAS, BL Lacertae had also 
been observed on a number of occasions, mostly with the full array. The data from those 
observations are also used in this work to establish a longer baseline. The total live exposure 
time (after quality selection) amounts to 20.3 hrs from 2010 September to 2011 November, 
with zenith angles ranging from 10 to 40 degrees. The source was not detected throughout 
the time period, except for one night on MJD 55740 (2011 June 28), when the automated 
realtime analysis revealed the presence of a rapidly flaring gamma-ray source in the direction 
of BL Lacertae. On that night, BL Lac ertae was observed only with three telescopes in the 



"wobble" mode (lAharonian et al.l 120011 ) with 0.5 degrees offset, because one telescope was 
temporarily out of commission. Starting at 10:22:24 UTC, two 20-minute runs were taken 
on the source under good weather conditions, with the zenith angle varying between 10 
and 13 degrees. No additional runs were possible due to imminent sunrise. The total live 
exposure time was 34.6 minutes. 
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The data were analyzed using the data analysis package described in ICogan (120081) . The 
analysis procedure includes raw data calibration, image par ameterization (lHillaslll985f ). event 
reconstruction, background rejection and signal extraction ( lDanielll2008f ). The standard data 
quality cuts (identical for the four- and three-telescope configuration), which were previously 
optimized for a simulated soft point source of ~6.6% of the Crab Nebula flux at 200 GeV 
and a photon index of 4, were applied to the shower images. The cuts used were: an 
integrated charge lower cut of 45 photoelectrons, a distance (between the image centroid 
and the center of the camera) upper cut of 1.43 degrees, a minimum number of pixels cut 
of 5 for each image, inclusive, mean scaled width and length cuts 0.05 < MSW < 1.15, 
and 0.05 < MSL < 1.3, respectively. A cut of 9 2 < 0.03 deg 2 on the size of the point- 
source search window was made, where 9 is the angle between the reconstructed gamma-ray 
direction and the direction to the source. A specific effective area corresponding to these 
cuts and the relevant array configuration was generated fro m simulations and was used to 



calculate the flux. The reflected-region background mo del (iBerge et al. 



for background estimation, a generalized method from iLi fc Mai (119831 ) was used for the 



20071 1 was applied 



calculation o f statistical s ignifi cance, and upper limits were calculated using the method 
described by iRolke et al.l (120051 ). The results were confir med by an ind ependent secondary 
analysis with a different analysis package, as described in IPanieJ ( 2008 ). 



2.2. High Energy Gamma Ray 

The Fermi Large Area Telescope (LAT) is a pair-conversion high-ener gy gamma-ray 



telesc ope covering an energy range from about 20 MeV to more than 300 GeV (lAtwood et al. 



20091 ). It has a large field-of-view of 2.4 sr, and an effective area of ~ 8000 cm 2 for > 1 GeV. 



In its nominal (survey) mode, the Fermi-LAT covers the full sky every 3 hours. 

During the time window when VERITAS detected a rapid flare on MJD 55740 (2011 
June 28), BL Lacertae was in the field of view of the LAT for about 16 minutes (MJD 
55740.431 - 55740.442). In analyzing the simultaneous LAT data, we selected Diffuse class 
photons with energy between 0.2 and 10 GeV in a 16° x 16° region of interest (ROI) centered 
at the location of BL Lacertae. Only events with rocking angle < 52° and zenith angle 
< 100° were selected. The data were processed using the publicly available Fermi-LAT tools 
(v9r23pl) with standard instrument response functions (P7S0URCE_V6). For such a short 
exposure, a very simple model containing the source of interest and the contribution of the 
galactic (using file gal_2yearp7v6_v0 . f its) and isotropic (using file iso_p7v6source.txt) 
diffuse emission was used. The contribution of the other known gamma-ray sources in the 
ROI is assumed to be negligible compared to that of BL Lacertae and the diffuse emission. 
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The model is fitted to the data using a binned likelihood analysis (gtlike), where the 
only free parameters are the spectral normalization and the power-law index of BL Lacertae. 
The contribution of the galactic and isotropic diffuse emission was fixed to a normalization of 
1.0, which is compatible with the values obtained when analyzing the same field of view dur- 
ing longer timespans. The results are used to construct an energy spectrum of BL Lacertae. 
We also performed an unbinned likelihood analysis and obtained similar spectral results. 

For comparison, we repeated the analyses for a longer period (of 24 hours) centered at 
the time of the VERITAS observations, as well as for times prior to the VERITAS-detected 
flare (between 2011 May 26 and 2011 June 26, or MJD 55707-55738). For the latter, we 
adopted a source model that incorporates all sources in the 2FGL catalog within the ROI 
and within 5 degrees of the ROI edges. The spectral results were extracted by adopting a 
custom spectral code (SED_scripts) available on the Fermi-LAT website. In all cases, the 
LAT spectrum of BL Lacertae can be well described by a power law, which justifies the 
assumption made in the likelihood analyses. 

A daily-binned light curve integrated above 0.1 GeV was derived covering the period 
MJD 55652-55949 (2011 April 01 - 2012 Jan 23) using the likelihood method described above. 
In each 1-day bin, the flux and the corresponding la error are calculated if the test statistic 
(TS) value is greater than 1, otherwise an upper limit is calculated. 



2.3. X-ray and Ultraviolet 

BL Lacertae was also ob served with the XRT and UVOT instruments on board the Swift 



satellite (IGehrels et al.ll2004j ) contemporaneously with the gamma-ray flare in 18 exposures 
between MJD 55704 (2011 May 23) and MJD 55768 (2011 July 26), including six ~2 ks 
Target of Opportunity (ToO) observations on six nights following the VHE flare on MJD 
55740. The combination of the X-ray telescope (XRT) and UV/optical telescope (UVOT) 
provided useful coverage in soft X-rays and UV, although none of the observations were 
simultaneous with the VERITAS observations during the flare. 

We analyzed the XRT data using the HEASOFT package (version 6.11). The event files 
are calibrated and cleaned using the calibration files from 2011 September 5. The data were 
taken in the photon-counting (PC) mode, and were selected from grades to 12 over the 
energy range 0.3-10 keV. Since the rates did not exceed 0.5 counts per second, pile-up effects 
were negligible. Source counts were extracted with a 20 pixel radius circle centered on the 
source, while background counts were extracted from a 40 pixel radius circle in a source-free 
region. Ancillary response files were generated using the xrtmkarf task, with corrections 
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applied for the point-spread function (PSF) losses and CCD defects. The corresponding 
response matrix from the XRT calibration files was applied. The spectrum was fitted with 
an absorbed power law model, allowing the neutral hydrogen (HI) column density (N H ) to 
vary. The best fitted value of Nh is (0.24 ± 0.01 ) x 10 22 cm" 2 , which is in agreement with the 
result of Nh 
val ue of Nh z 



0.25 x 10 22 cm 2 presented by Ravasio et al. ( 20031 ). but is larger than the 
0.18 x 10 22 cm -2 from the Leiden/Argentine/Bonn (LAB) survey of galactic 



HI (IKalberla et all 120051 ). 



The UVOT cycled through each of the optical and the UV pass bands V, B, U, UVW1, 
UVM2 and UVW2. Data were taken in the image mode discarding the photon timing 
information. Only data from UVW2 band are shown in this work; the other bands roughly 
track UVW2. The photometry w as c omputed using a n ape rture of 5" following the general 
prescription of iPoole et al.l (120081 ) and lBreeveld et al.l (l2010l ). Contamination by background 
light arising from nearby sources was removed by introducing ad hoc exclusion regions. 
Adopting the N H val ue provided by the XRT a nalysis and assumin g; E(B — V) = 0.34 mag 
(IMaesano et al.lll997l ). we estimated Ry = 3.2 (IGiiver fc Ozelll2009l ). Then, the optical/UV 
galactic extinction coefficients were applied ( Fitzpatrickl 1999). The host g alaxy contribution 



has been estimated using the PEGASE-HR code (ILe Borgne et al.l 120041 ) extended for the 
ultraviolet UVOT filters. Moreover, there is no pixel saturation in the source region and no 
significant photon loss. Therefore, it is possible to constrain the systematics to below 10%. 



2.4. Optical 



As part of the Steward Observatory spectropolarimetric monitoring project (ISmith et al. 



20091 ). BL Lacertae was observed regularly with the 2.3m Bok Telescope and the 1.54m 
Kuiper Telescope in Arizona. Measurements of the V-band flux density and optical linear po- 
larization are from the Steward Observatory public data archive Qhttp : // j ames . as . arizona . edu/~ psmith 



The d ata were reduced and calibrated following the procedures described by ISmith et al. 
(120091 ). We note that there is a 180-degree degeneracy in polarization angle, so we shifted 



some polarization angles by 180 degrees to minimize the change between two consecutive 
measurements. No corrections to the data have been made for the contribution from the 
host galaxy, or interstellar polarization, extinction and reddening. However, these issues 
have little effect on variability studies. 
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2.5. Radio 



BL Lacertae was observed with the VLBA at 43 GHz, roughly once a month, as part of 



the monitoring program of gamma-ray bright blazars at Boston University (http : / /www . bu . edu/blazars/^ 
Two extra epochs of imaging were added via Director's Discretionary Time on 2011 July 6 
and 29. The data were correlated at the National Radio Astronomy Observatory in So- 
corro, NM, a n d the n analyzed at Boston University following the procedures outlined by 
Jorstad et al.l (120051 ). The calibrated total and polarized intensity images were used to 
investigate the jet kinematics and to calculate the polarization parameters (degree of polar- 
ization p and position angle of polarization x) f° r the whole source imaged at 43 GHz with 
the VLBA and for individual jet components. The uncertainties of polarization parameters 
were computed based on the noise level of total and polarized intensity images and do not 
exceed 0.6% and 3.5 degrees for degree of polarization and position angle of polarization, 
respectively. 

BL Lacertae is also in the sample of the Monitoring Of Jets in Active galactic nuclei 
with VLBA Experiments (MOJAVE) program. For this work, we only used results from 
polarization meas urements at 15.4 GHz. The data reduction procedures are described by 
Lister et al.l ( 120091 ). Briefly, the flux density of the core component is derived from a Gaussian 
model fit to the interferometric visibility data. Polarization properties of the core are then 
derived by taking the mean Stokes Q and U flux densities of the nine contiguous pixels 
that are centered at the Gaussian peak pixel position of the core fit. The results include 
fractional linear polarization, electric vector position angle (note the 180-degree degeneracy), 
and polarized flux densities. The flux density has an uncertainty of ~ 5%, and the position 
angle of polarization has an uncertainty of ~ 3 degrees. 

For better sampling, we used data from blazar monitoring programs with the Owens 
Valley Radio Observatory (OVRO) at 15.4 GHz, with the Metsahovi Radio Observatory 
(MRO) at 37 GHz, and with the Submillimeter Array (SMA) at 230 and 350 GHz, respec- 
tively. The OVRO 40 m uses off-axis dual-beam optics and a cryogenic high electron mobility 
transistor (HEMT) low-noise amplifier with a 15.0 GHz center frequency and 3 GHz band- 
width. The two sky beams are Dicke-switched using the off-source beam as a reference, and 
the source is alternated between the two beams in an ON-ON fashion to remove atmospheric 
and ground contamination. Calibration is achieved using a temperature-stable diode noise 
source to remove receive r gain drifts and th e flux density scale is derived from observations 
of 3C 286 assuming the iBaars et al.l (119771 ) value of 3.44 Jy at 15.0 GHz. The systematic 
uncertainty of about 5% in the flux density scale is not include d in the error ba r s. Co mplete 
details of the reduction and calibration procedure are found in iRichards et al.l (120 111 ). 



The 37 GHz observations were made with the 13.7 m diameter Metsahovi radio telescope, 
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which is a radome-enclosed paraboloid antenna situated in Finland (24 23' 38" E, +60 13' 
05"). The measurements were made with a 1 GHz-band dual beam receiver centered at 
36.8 GHz. The observations are ON-ON observations, alternating the source and the sky in 
each feed horn. A typical integration time to obtain one flux density data point is between 
1200 s and 1400 s. The detection limit of the telescope at 37 GHz is on the order of 0.2 
Jy under optimal conditions. Data points with a signal-to-noise ratio < 4 are treated as 
non-detections. The flux density scale is set by observations of DR 21. Sources NGC 7027, 
3C 274 and 3C 84 are used as se condary calibrato r s. A detailed description of the data 



reduction and analysis is given in iTerasranta et all ( 119981 ). The error estimate in the flux 



density includes the contribution from the measurement rms and the uncertainty of the 
absolute calibration. 

Observations of BL Lacertae at frequencies near 230 and 350 GHz are from the Submil- 
limeter Array (SMA), a radio interferometer consisting of eight 6-m diameter radio telescopes 
located just below the summit of Mauna Kea, Hawaii. These data were obt ained and call 



brated as part of the normal monitoring program initiated by the SMA (see iGurwell et al. 



20071 ). Generally, the signal-to- noise ratio of these observations exceeds 50 and is often well 
over 100, and the true error on the measured flux density is limited by systematic rather 
than signal-to-noise effects. Visibility amplitudes are calibrated by referencing to standard 
sources of well-understood brightness, typically solar system objects such as Uranus, Nep- 
tune, Titan, Ganymede, or Callisto. Models of the brightness of these objects are accurate 
to within around 5% at these frequencies. Moreover, the SMA usually processes only a single 
polarization at one time, and there is evidence that BL Lacertae in 2011 exhibited a fairly 
strong (~15%) linear polarization. For a long observation covering a significant range of 
parallactic angle, the effect of the linear polarization would be largely washed out, providing 
a good measure of the flux density. However, not all observations of BL Lacertae covered 
a significant range of parallactic angle, and thus in some cases we would expect a potential 
absolute systematic error up to 10%. In most cases, we expect that the total systematic 
error is around 7.5%. 



3. Results 

3.1. Gamma Ray Properties 

The VERITAS analysis showed an excess of 212 7-like events, corresponding to 11.0 ± 
0.8 7/min and a 21. la detection of BL Lacertae in the first observation run on MJD 55740 
(2011 June 28), with an effective exposure of 19.3 minutes starting at 10:22:24 UTC. The 
second run, with an effective exposure of 15.3 minutes, yielded an excess of only 33 7-like 
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events, corresponding to a 4.1cx detection. The VERITAS analysis of 19. 7- hour data from 
2010 September to 2011 November, excluding the two flaring runs, showed an excess of 21 
7-like events, and a statistical significance of 0.28o\ 

Focusing on the two flaring runs, we produced a light curve with 4-minute bins as shown 
in the inset of Fig. [TJ The fluxes were computed with a lower energy threshold of 200 GeV. 
The observations missed the rising phase of the flare. In 4-minute bins, the highest flux that 
was measured is (3.4 ± 0.6) x 10~ 6 photons m _2 s _1 , which corresponds to about 125% of 
the Crab Nebula flux above 200 GeV, as measured with VERITAS. To quantify the decay 
time, the light curve was fitted with an exponential function I(t) = Iq x exp(— t/ra), and 
the best-fit decay time was = 13 ± 4 minutes. 

To determine the position of the flaring gamma-ray source, we fitted a 2-d Gaussian 
function to the uncorrelated map (binned to 0.05°) of excess events (after acceptance correc- 
tion) from both runs. The best-fit right ascension and declination ( J2000) are 22 /l 02 m 37 s and 
+42°15'25 ", respectively, with a statistical uncertainty of ~0.01° along both directions. The 
source is thus named as VER J2202+422. According to the Simbad database, BL Lacertae 
is the only object within a radius of 2'. 

Using data from the first flaring run, we extracted a gamma-ray spectrum (Fig. [2]). It 
can be fitted with a power law, 

dN/dE = (0.58 ± 0.07) x 10~ 9 x ( J E/0.3TeV) ( - 3 ' 6±0 - 4) cm- 2 S - 1 TeV~ 1 . (1) 

Also shown in the figure is the gamma-ray spe ctrum of BL Lacertae obtained with MAGIC 



during a lower- flux state (lAlbert et al.l 120071 ). The two gamma-ray spectra have nearly 



the same sl ope. This is in cont rast with the typical spectral hardening trend of a flaring 



blazar (e.g.. iGiommi et al.lll990l ). It may reflect the fact that in both cases the TeV gamma 
rays fall on the steeply falling part of the high-energy SED peak, which might not be sensitive 
to flux changes. A spectrum was also constructed from both runs and fitted with a power 
law: 

dN/dE = (0.30 ±0.03) x 10~ 9 x (E/0.3TeV) ( - 3 ' 8±0 - 3) cm- 2 s- 1 TeV- 1 . 

To better constrain the gamma-ray SED, we plotted the Fermi-LAT spectra of the source 
averaged over several time periods (16 minutes, one day and one month) along with the 
VERITAS spectra, in Fig. |3j The VERITAS spec tra during the flare both with and without 



extragalactic background light (EBL) corrections (IDommguez et al.ll201ll ) are shown, as well 
as the 95% confidence upper limits from 14 observations in one month before the flare. The 
best-fit LAT results for the 16 minutes simultaneous with VERITAS converge to -Fo.2-iogcV — 
(2.1 ± 0.9) x 10 _6 cm~ 2 s _1 with a spectral index of T = 1.6 ± 0.4, with a test statistic of 35. 
The 16-minute exposure is very short for the LAT, so the uncertainty is large. Together, the 
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simultaneous VERITAS and LAT spectra show that the gamma-ray SED peak probably lies 
between 10 and 100 GeV. Moreover, the LAT results provide evidence for spectral hardening 
during the VERITAS flare, with the best-fit photon index changing from about 2.12 ± 0.05 
to 1.6 ± 0.4 in the LAT band; note, however, the large uncertainties. 



3.2. Multiwavelength Properties 

Fig. H] shows the multiwavelength light curves of BL Lacertae. T he TeV gamma-ra y 



flare occurred when the source was active and variable at GeV energies (jCutini et al.ll201ll ). 
Although the LAT light curve shows variability on a timescale of days, no rapid flaring on 
shorter timescales is apparent. On the other hand, the LAT could have missed a flare as 
rapid as the TeV flare, due to the lack of statistics. There was no apparent activity in the 
soft X-ray either, but changes in flux may have occurred at UV and optical wavelengths 
(although there was no UVOT coverage during the TeV gamma-ray flare). The XRT spec- 
trum was relatively hard, with a photon index of ~1. 8, compared to th at obtained during 



the Fermi-L AT campaign during a low state in 2008 (lAbdo et al.ll201ll ). There was no ap- 
parent variation in the radio flux from the source at the time of the TeV gamma-ray flare. 
It is interesting to note that the LAT caught a very intense rapid flare earlier (around MJD 
55710) that was accompanied by similar flaring activities at soft X-ray, UV, and optical 
wavelengths. Unfortunately, the source was not observed with VERITAS at that time. 

Extending the Fermi-LAT, optical, and radio light curves to later times, we clearly see 
an intense flare that occurred at 15.4 GHz, 37 GHz, and 230 GHz, as shown in Fig. [5j about 
four months after the TeV gamma-ray flare. Although the elevated flux is also evident at 
350 GHz, the flare is poorly sampled. The well-sampled Fermi-LAT light curve indicates 
some elevated and variable GeV emission in 2011 November. However, the presence of 
similar GeV variabilities from 2011 May to the end of the year makes it difficult to establish 
a correlation between GeV and radio bands. We cross-correlated the light curves at the four 



radio frequencies, using the z-transformed discrete correlation function (ZDCF) ([Alexander 



19971 ). The results are shown in Fig. |6l indicating high degree of correlation among the 
bands. From the ZDCFs, the corresponding time lags were measured, using a publicly 
available likelihood code (PLIKE), and are plotted against v~ l in Fig. [7J 

Fig. |H] shows results from radio polarization measurements. Although there is no sig- 
nificant variation in the average polarization fraction, the average polarization angle of the 
core appears to change before and after the TeV gamma-ray flare. However, the polariza- 
tion angles for VLBA 15.4 GHz and 43 GHz do not agree with each other in earlier epochs 
(before the TeV gamma-ray flare). This discrepancy is likely due to the combination of the 
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emergence of a new component, the Faraday rotation and the difference in beam size at 
the two frequencies. At the core, the Fa raday rotation can be significant for BL Lacertae 



( IGabuzda et all 12006c ; iJorstad et al.l 120071 ). mostly affecting the 15.4 GHz measurements. It 



is also worth noting that the effects can be variable on timescales of months. 

The emergence of a new component is strongly supported by the VLBA observations at 
43 GHz. Fig. [9] shows a series of high-resolution images of BL Lacertae around the time of 
the gamma-ray flare, with the measured polarization flux and angle indicated. The results on 
the long-term VLBA monitoring observations will be presented elsewhere ( Jorsted, Marscher 
et al., in prep). The new knot, Kll, is discerned from the core in the images by its different 
polarization position angle x (20° for Kll compared with 44° for the core), even before it is 
clearly seen in the total intensity contours in the 2011 July 29 image. The proper motion of 
Kll is 0.72 mas yr^ 1 , which corresponds to an apparent speed of 3.6c, with uncertainties of 
order 20% owing to the short time interval over which the trajectory is followed. Although 
they have lower resolution, the MO JAVE images, as shown in Fig. [101 also indicate a change 
in the polarization of the core before and after the TeV gamma-ray flare but not in the 
downstream jet polarization, lending further support for the emergence of a new component 
associated with the gamma-ray flare. 

Fig. [8] also shows results from optical polarization measurements. Again, the polarized 
flux does not vary significantly before and after the gamma-ray flares. However, changes in 
optical polarization angle are significant around the times of both GeV and TeV gamma-ray 
flares and between. Over the 4-day period that included the VERITAS flare, the optical 
polarization position angle changed by a minimum of 38.8 deg (between MJD 55738 and 
55739), -31.2 deg (between MJD 55739 and 55740), and 88.8 deg (between MJD 55740 and 
55741). Therefore, at a minimum, the optical polarization angle was changing by more than 
one degree per hour. A similar pattern is seen for the Fermi LAT flare earlier. 



4. Discussion 

For the first time, a rapid (minute-scale) TeV gamma-ray flare is seen from BL Lacertae - 
this is the first such flare from an LBL. It fills an important gap between similar phenomenon 
observed in FSRQs and HBLs. Independent of any models, the measured decay time of the 
flare (r^) requires that the size of the emitting region must be very small, 

R < CT d 5/{\ + z)^ 2.2 x 10 13 5 cm, (2) 

where z is the redshift of the source (z = 0.069) and 5 is the Dopper factor of the jet, 

5= [T(l - PcosO)}- 1 ; 
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T is the bulk Lorenz factor of the jet; and 6 is the angle between the axis of the jet and the 
line of sight. 

Another constraint for the Doppler factor is derived from the fact that the emitted 
gamma r ays must escape 77 pair p roduction in the emitting regions. The optical depth is 
given by ( iDondi fc Ghisellinilll995l ) 



T- 



\2a 
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5 



-(4+2o) 

5hc 2 T1/2 



(3) 



where is the luminosity distance; T1/2 is the doubling (not e-folding) time of the flare; the 
fiducial frequency is defined as 



v = - 

V 



m e c 



where v is the observed frequency of a gamma-ray photon; F(u ) is the observed flux at 
is . We should note that the frequency of the target photon is fo[<5/(l + z )] m the frame 
and ^o[^/(i + z)} 2 in the observer frame. Eq. [3] is derived under the assumption that the 
energy spectrum of the source is approximated by introducing a power law between uq and 
u [5/(l + z)} 2 , with a spectral index of a. The fact that we detect TeV gamma-ray emission 
implies that the optical depth r 77 cannot be too large. Requiring r 77 < 1 leads to 

1/(4+2q) 



5 > 



a T d 2 L (l + z) 2a F(u ) 



5hc 2 



Tl/2 



(4) 



The calculation assumes that gamma rays and target photons are both isotropic in the jet 
frame, so is, strictly speaking, only applicable if the gamma rays are produced via synchrotron 
self-Compton (SSC) scattering. For BL Lacertae, SSC may be a good approximation, as 
BLRs are quite weak. In this case, di ~ 311Mpc and T1/2 ~ 9 minutes, assuming Q m = 0.27, 
Q A = 0.73, and H = 70 km s^Mpc" 1 Larson et al.ll201lh . At hv w 0.9 TeV, which is about 
the highest energy of all gamma rays detected within the source region, we have i/ ~ 7x 10 13 
Hz. Unfortunately, we did not have simu l taneo us IR coverage during the gamma-ray flare. 



Interestingly, according to iRaiteri et al.l ( 120091 ). the IR flux of BL Lacertae did not vary 



significantly (within a factor 2) durin g; their long-term monitoring (for over 150 days) in 



2007-2008 (see also lAbdo et al.ll201ll ). It is also worth noting that the synchrotron SED 



peak of BL Lacertae lies in the near-IR band, and the archival SEDs between near-IR and 



X-ray can be rou ghly described by a power law (e.g.. lBottcher et al.ll2003l ; IRaiteri et al.ll2010 



Abdo et al.ll201ll ). The spectral index (a) varied in the range 1.34-1.40 in the fre quency range 



7 x 1 13 -10 17 Hz. Takin g F(vq) and a from the archival SEDs of BL Lacertae (IRaiteri et al. 



20091 ; lAbdo et al.1 1201 ll ). from Eq. H] we found that the lower limit on 5 lies in the range 
13-17. 



- 15 - 



The derived lower limi ts on 5 can be comp ared with other estimates. From the radio 
variability of BL Lacertae, iHovatta et al.l (120091 ) derived a value of 5 = 7.3. However, the 
uncertainty is expected to be large due to a number of assumptions involved in the analysis 
(especially in rela tion to the int r insic brightness temperature). On the other hand, using a 
different method, iJorstad et al.l (120051) arrived at a v alue o f 5 = 7.2 ±1.1 for different jet 



components, in good agreement with IHovatta et al.l (120091 ) . These values are significantly 
below the lower limits imposed by gamma-ray observations, perhaps implying differences 
between radio and gamma-ray emitting regions in the jet or a gamma-ray optical depth of 
r 77 ~ 40. It remains to be seen whether such a strong attenuation of TeV gamma rays can 
be accommodated theoretically. The efforts to model the broadband SEP o f BL Lacertae 
have generally led to Do ppler factors larger than 7 (e.g., lAbdo et al.l 1201 ll ; see, however, 
Bottcher fc Reimerll20oi ). 



Rapid TeV gamma-ray flaring was first observed in HBLs. It was recognized immedi- 
ately that the requisite (large) Doppler factor would be p roblematic, becau se no superlu- 
minal motion had ever been seen in any of these sources ( iPiner et al.l 120081 1 . This led to 
the suggestion of a s tratified structure of the jet that consists of a fast-moving spine and 
slow-moving sheath ( iGhisellini et al.l 120051 1 . The high-resolution polarization maps of the 



TeV gamma-ray HBLs have provided some evidence for such a configuration ( IPiner et al. 



20081 ). However, for BL Lacertae, the polarization measurements do not show any stratifica- 
tion of the jet (for example, see Fig. [9]). There is certainly no evidence for a slowly-moving 
sheath. Alternatively, the large Doppler factor might imply that gamma-ray production 
occurs in a region upstream of what is observed with the VLBA. Deceleration could explain 
the discrepancy between the values of the Doppler factor inferred from gamma-ray and ra - 
dio observations (IGeorganopoulos fc Kazanasll2003l ; lLevinsonll2007l : IStern &: Poutanenll2008l ). 
However, pushing the gamma-ray production region too close to the central black hole would 
be problematic for BL Lacertae and, even more so, for PKS 1222+21, as attenuation due to 
radiation from the BLRs would be strong. These difficulties might be alleviated in models 
that invoke subregions inside the jets that are fast moy i ng and also sufficiently far from the 
black hole (iGiannios et al.ll2009l ; lNarayan fc Piranll2012t iNalewajko et al.ll2012t see, however, 
Tavecchio et al.l 1201 ll ). 



Of particular significance is that our multiwavelength observations of BL Lacertae might 
link the emergence of a superluminal knot (Kll in Fig. [9]) with the TeV gamma-ray flare. 
The former is directly seen in the VLBA 43 GHz images, although there is a large gap in the 
coverage around the time of the TeV gamma-ray flare. The VLBA 15 GHz observations also 
show changes in the polarization angle, which supports the emergence of a new component. 



Based on earlier VLBA imaging, iMarscher et al. 



(120081 ) argued that the core is a standing 



shock located well downstream (by ~ 1 pc) of the black hole. Their model also describes a 
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helical magnetic field configuration upstream of the radio core, which the radiating plasma 
follows. This is now supported by the observed pattern of change in the optical polarization 
that coincides with the TeV gamma-ray flare. The new superluminal knot seems to have 
passed through the core on MJD 55711±15 (2011 May 30, when the brightness centroids 
of the knot and core coincided), close to the time when a rapid flare was seen with the 
Fermi-LAT, Swift XRT and UVOT, and the Steward Observatory. 



In the model of iMarscherl (120121 ). the radio core is a conical shock that ends in a small 
shock oriented transverse to the jet axis (a Mach disk). The slow but highly compressed 
plasma in the Mach disk provides a highly variable local source of seed photons for inverse- 
Compton scattering by electrons in the faster plasma that passes across the conical shock. 
If a region of especially high density of relativistic electrons passes through the core, it 
can cause a sharp flare at gamma-ray energies and appear as a superluminal knot at radio 
frequencies. Although the angular resolution of the VLBA is insufficient to measure the 
angular size of the knot during the observations, it is likely to have a diameter ~ 0.07 mas 
assuming that its brightness temperature is close to the value of ~ 5 x 10 10 K needed for 
equip artition between the energy density in relativistic electrons and that in the magnetic 
field (jReadheadlll994l ) . In this case, the knot interacted with the core over a period of 70 ± 15 
days centered on MJD 55711 (i.e., from late April 2011 until early July 2011). Therefore, 
the knot would be near the end of the core region when the TeV gamma-ray flare erupted. 

Alternatively, the burst of TeV gamma rays may be produced far upstream of the radio 
core (closer to the supermassive black hole), related to the emergence of a high-density re- 
gion. As the reg ion moves downstream, and along the helical magnetic field (as postulated by 



Marscherl (120121 )). it produces polarized optical emission with a characteristic variation pat- 
tern of the polarization angle. When it becomes optically thin to synchrotron self-absorption, 
still further downstream, it is seen at successively longer wavelengths. The observed radio 
flare may be a manifestation of it. A delay of the radio flare by four months, with respect 
to the gamma-ray precursor, is in line with the fact th 
erally lags the gamma-ray variability by 1-8 months 



months (e.g., 


Lahteenmaki & Valtaoia 


2003; 


Leon-Tavares et al. 


2012 


). Theoretically, the 



optical-depth effect should lead to a v 1 dependence of the time lag, as the core remains 
optically thick to synchrotr on self-absorption up to a distance (from the black hole) r c oc z/ _1 



(IBlandford &: Konigllll979l ). Unfortunately, the measurements (as shown in Fig. [7]) are not 



of sufficient quality to confirm such a frequency dependence. 

The lack of similarly rapid change of significant amplitude at other wavelengths is likely 
due to inadequate sampling. In other words, the TeV gamma-ray flare is so rapid that 
pointed instruments were unlikely to be observing the source at the right time, while for 
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other instruments (e.g., the Fermi LAT) it is difficult to accumulate adequate statistics. 
Nevertheless, around the time of the TeV gamma-ray flare, there is evidence for flux varia- 
tions at optical and UV wavelengths, which would represent a response of the synchrotron 
emission to the VHE gamma-ray flaring. 
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Fig. 1. — TeV gamma-ray light curve of BL Lacertae (> 200 GeV). When the source was 
not significantly detected, 99% confidence upper limits are shown. The upper limits were 
derived by combining data from all observation runs for each night, but for the night of the 
flare, the fluxes derived from the two individual runs are shown separately. The inset shows 
the flare in detail, in 4- minute bins for the first run, and one 16-minute bin for the second 
run, with minute indicating the start of the first run. The dashed line shows the best fit 
to the profile with an exponential function (see text). 
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Fig. 2. — TeV gamma-ray spectrum of BL Lacertae. The VERITAS data points are shown 
as red open diamonds, along with the best-fit power law (solid line). For comparison, we 
also show the published MAGIC spectrum of the source as blue open squares, along with 
the best-fit power law (dashed line). The two power laws have comparable slopes. 
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Fig. 3. — Broadband gamma-ray spectrum of BL Lacertae. For VERITAS, the flare spectrum 
with EBL correction is shown as blue points with uncertainties, and the spectrum without 
EBL correction is shown as lighter blue points. Also plotted are the 95% confidence upper 
limits that were derived from 14 VERITAS observation runs in the month prior to the flare. 
For the Fermi-LAT data, the simultaneous, 1-day, and pre-flare spectra are shown in blue, 
light blue, and gray. The pre-flare spectrum was derived from LAT observations taken over 
a period of one month before the TeV gamma-ray flare. The shaded areas show the la 
confidence interval of the overall Fermi derived spectra, independent from the data points. 
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Fig. 4. — Multiwavelength light curves of BL Lacertae. Fermi-LAT la upper limits are 
shown as arrows. The dotted line indicates the time of the TeV gamma-ray flare seen with 
VERITAS, while the dashed line shows the time of a rapid GeV gamma-ray flare seen with 
the Fermi-LAT. 
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Fig. 5. — Extended Fermi-LAT, optical, and radio light curves of BL Lacertae. As in Fig. HJ 
the dotted line indicates the time of the VERITAS flare, and the dashed line shows the time 
of the Fermi-LAT flare. 
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Fig. 6. — Cross-correlation functions derived from the light curves of BL Lacertae: top 
between the SMA 350 GHz and OVRO 15 GHz bands; middle between the SMA 230 GHz 
and OVRO 15 GHz bands; and bottom between the Metsahovi 37 GHz and OVRO 15GHz 
bands. Positive delays indicate "lead" with respect to the reference (OVRO 15 GHz) band. 
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Fig. 7. — Time delays of the radio flare from BL Lacertae. The time delays with respect 
to the OVRO band were determined from a likelihood code PLIKE. Positive delays indicate 
"lead" with respect to the OVRO band. The dashed line is drawn to guide the eye. See the 
text for discussion. 
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Fig. 8. — Polarized optical and radio emission from BL Lacertae. VERITAS and Fermi- 
LAT light curves are also shown for comparison. The optical measurements were made 
at the Steward Observatory, while the radio measurements were made with the VLBA at 
15.4 GHz (black dots) and 43 GHz (red triangles). The radio electric vector position angle 
has an uncertainty of about ±3 degrees, and the polarization fraction has an uncertainty of 
about 5%. The dotted line indicates the time of the VERITAS flare, and the dashed line 
shows the time of the Fermi-LAT flare. 
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Fig. 9. — 43 GHz VLBA images of BL Lacertae at four epochs around the time of the TeV 
gamma-ray flare. The images are convolved with a circular Gaussian function (represented 
by the circle in the bottom- left corner) that has a full width at half maximum of 0.1 mas (i.e., 
~0.15 pc at the distance of 311 Mpc), the approximate resolution of the longest baselines of 
the array. Contours correspond to total intensity, with levels in factors of 2 from 0.25%, plus 
an extra contour at 96%, of the peak intensity of 2.16 Jy beam -1 . Color represents linearly 
polarized intensity, with maximum (black) of 0.103 Jy beam -1 followed by red, blue, yellow, 
and white (no polarization detected). Red lines mark the position of the assumed stationary 
core and the superluminally moving knot Kll, each of which has a distinct polarization 
position angle. 
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Fig. 10.— MOJAVE 15.4 GHz VLBA images of BL Lacertae at three epochs in 2011, 
showing a change in core polarization after the 2011 June 28 TeV flare. The images on the 
left show total intensity contours, with electric polarization vectors overlaid in blue. The 
images on the right show total intensity contours with fractional linear polarization in color. 
The polarization color scale ranges from to 50%. The images have been convolved with the 
same Gaussian restoring beam having dimensions 0.89 mas x 0.56 mas and position angle 
—8 degrees. The base contour levels in each image are 1.3 mJy beam -1 in total intensity 
and 1 mJy beam -1 in polarization. The angular scale of the image is 1.29 pc mas -1 . 



